Spin polarization and edge currents in topological superconductors 
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The spin character of topologically protected edge states in noncentrosymmetric superconductors is investi- 
gated through exact diagonalization of Bogoliubov-de Gennes Hamiltonians. We find that due to the spin-orbit 
coupling and parity-mixing, the edge states exhibit an intricate spin structure, which depends on the edge mo- 
mentum. If the superconductor is brought into contact with a ferromagnetic insulator, the exchange field cou- 
pling to the spin-polarized edge states gives rise to a spontaneous edge current. We study the variation of the 
zero-temperature edge current on the magnetization direction and the singlet-triplet mixing in the superconduc- 
tor, and show that it can exhibit singular dependence on the exchange-field strength. 
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Introduction. — The realization that effectively noninteract- 
ing insulators and superconductors can be systematically clas- 
sified in terms of their topological properties has inspired an 
intense research effort EHD . An important class of topologi- 
cal systems are noncentrosymmetric superconductors (NCSs), 
characterized by strong antisymmetric spin-orbit coupling 
(SOC) and a mixing of spin-singlet and spin-triplet pair- 
ing 0. Remarkably, both fully gapped and nodal NCS can 
have nontrivial topology, which is manifested by topologically 
protected edge states via the bulk-boundary correspondence. 
In a fully gapped NCS, the edge states are dispersing Majo- 
rana modes (UEirEI, whereas a nodal NCS exhibits nondis- 
persing zero-energy edge states (flat bands) [ 9 - 14 1. There 
are several candidate materials for topological NCSs, such as 
CePt 3 Si E2, Li 2 Pt 3 B ED, and BiPd G3- Furthermore, it 
has been proposed to engineer an NCS by exploiting the prox- 
imity effect in a material with strong SOC placed in contact 
with a topologically trivial superconductor lfT8l - [2TI , or in the 
metallic layer realized at oxide interfaces l22ll23ll . 

In this Letter, we examine two contrasting examples of 
NCSs: (i) a fully gapped system with helical Majorana edge 
states and (ii) a nodal state with nondegenerate flat bands. We 
demonstrate that these edge states are strongly spin polarized, 
which originates from both the SOC and parity mixing in the 
superconductor. If the superconductor is brought into contact 
with a ferromagnetic insulator, the coupling to the exchange 
field alters the dispersion of these subgap states so that they 
carry a charge current along the edge. Besides the edge- state 
contribution, there is also a contribution to the edge current 
from bulk continuum states. We find that the properties of the 
edge current reflect the topological character of the NCS: For 
example, the topologically protected flat bands of the nodal 
NCS are responsible for a strong edge current with singular 
dependence on the exchange-field strength. 

Model Hamiltonian. — Quasiparticle motion in an NCS is 
described by the Bogoliubov-de Gennes Hamiltonian H = 
E k ^Ic^k^k, with ^ k = (ckt, 014 > c -kf c -k|) T and 

_ fe^a - Al k • cr A k \ 

Hk -{ At _, kao + Al k . ( r*J' (1) 



where lk = x sin k y — y sin k x , = t (cos k x + cos k y ) — /i, 
and Ak = (^k^o + dk • cr) (ia y ). Here, c\. a denotes the 
electron creation operator with spin cr and momentum k, 
Ii represents the chemical potential, t is twice the nearest- 
neighbor hopping integral, A stands for the SOC strength, 
and cr denotes the vector of Pauli matrices. The even-parity 
spin-singlet and odd-parity spin-triplet components of the su- 
perconducting order parameter are parametrized as ^k = 
A /(k) q and d k = A /(k) l k (1 - q), respectively. With 
the parameter q, the NCS can be tuned from purely spin- 
triplet (q = 0) to purely spin- singlet (q = 1) pairing. The 
structure factor f(k) allows to study different orbital- angular- 
momentum pairing states. In the following we consider two 
cases: (i) /(k) = 1 for an NCS with (s+p)-wave pairing sym- 
metry and (ii) /(k) = sin k x sin k y for the so-called (d xy +p)- 
wave pairing state. In our numerical calculations we choose 
(t,/z,A,A ) = (4.0,4.0,-2.0,0.5). 

Edge states. — To examine the different topological proper- 
ties of (s+p)-wave and (d xy +p)-wave pairing states, we study 
first the subgap states at the edge of the NCS, both with 
and without a proximity-induced exchange field H ex acting 
on the leading edge. To that end, we compute the spin- 
and momentum-resolved density of states (DOS) of Hamilto- 
nian ([T]) in a ribbon geometry with edges perpendicular to the 
(10) direction. The momentum-resolved DOS and the spin- 
resolved DOS in the n-th layer from the edge are given by 

p n (E, k y ) = -i- Im [Tr {G™(n> n; E)}] , (2a) 

ME, k y ) = - A im [ Tr G^\n, n; E)}] , (2b) 

respectively, where = diag (a M , — [a^]*). Here, 

G^® = (E + if] — i^ 0) ) _1 denotes the zero-temperature 

lattice Green's function of the ribbon with = 

± J dk x e^^-- / )i7k + ^^l^Adiag(H ex • ^ 7,-[H ex • ^ 7] T ). 
The expressions in Eqs. ^ are evaluated by taking a ribbon of 
width up to N = 10 3 and an intrinsic broadening rj = 0.005. 

We first discuss the edge DOS for H ex = 0. Figure [T] shows 
the local DOS at the (10) edge integrated over the ten outer- 
most layers, i.e., Y^f=i Pn{E, k y ), for the triplet-dominated 
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FIG. 1. (a) and (b): Momentum-resolved DOS on a log scale for 
the ten outermost layers at the (10) edge of an (s+p)- wave and a 
(d xy +p)-W2Lve NCS, respectively, (c), (d) x component and (e), (f) z 
component of the spin-resolved DOS for the outermost layer at the 
(10) edge in units of ft/10 on a linear scale. In all cases we choose 
q = 0.25, representative for the triplet-dominated pairing state. 



pairing state. Ten layers are on the order of the superconduct- 
ing coherence length, i.e., the distance over which the edge 
states are localized. Here, we present the edge DOS only 
for the triplet-dominated phase; the edge DOS for the singlet- 
dominated phase is included as Fig.[6]in the Supplemental Ma- 
terial (SM) for completeness. Figures [T] and [6] reveal both the 
subgap states (dark red) and the bulk continuum states (orange 
and yellow). In the parameter ranges q < q c ^ ~ 0.472 and 
q > q c jj ~ 0.583, the (s+p)- wave state is a fully gapped su- 
perconductor in symmetry class DHL For q < q c ^, the super- 
conductor has a topologically nontrivial character with a non- 
zero Z 2 topological number (2l|8]|25]]. By the bulk-boundary 
correspondence, we accordingly find helical Majorana subgap 
states, see Fig.[TJa) (6H3- For q > q Cy u, the system is topo- 
logically trivial and the subgap states hence disappear l24l . 
The (d xy +p)-wave NCS, in contrast, is nodal and therefore 
does not have a quantized global topological number for any 
q. However, by treating every point in the edge Brillouin zone 
(BZ) as the edge of a one-dimensional system, one can define 
a momentum-dependent winding number W(\o)(k y ), which 
only changes across projected nodes of the bulk gap fT0lfT3ll . 
Provided that q £ [q c ,L,q c ,u], one finds W(\o)(k y ) — ±1 for 
\k y \ e [k Fi +,k F -] ~ [0.3527T, 0.648tt]. That is, when k y 
lies between the projected edges kp,+ and hp,- of the two 
spin-orbit- split Fermi surfaces, the winding number W(\o)(k y ) 
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FIG. 2. (Color online) Surface DOS on a log scale for the ten 
outermost layers at the (10) edge of an (s+p)- wave and a (d xy +p)- 
wave NCS in the presence of an exchange field in the leading edge 
with magnitude h ex = 0.4. (a) and (c): exchange field along the 
x axis, (b) and (d): exchange field along the z axis. As in Fig. [T] 
we choose q = 0.25, representative for the triplet-dominated pairing 
state. 



takes on the nonzero values ±1, which ensures the existence 
of nondegenerate zero-energy flat bands at these momenta 
k y (9). For \k y \ < fcp,+> ° n the other hand, there are topolog- 
ically trivial dispersing states for q < q c ^, and doubly degen- 
erate zero-energy states with W^(k y ) = ±2 for q > q c jj, 
see Figs.[TJb) and|6jb). 

In the lower two rows of Fig. [T] we present the spin- 
resolved edge DOS p% =1 (E,k y ) in the absence of an ex- 
change field. The corresponding plots for the singlet-domi- 
nated phase are given in Fig. [6] in the SM. Interestingly, both 
the continuum and subgap states show a strong polarization 
in the xz spin-plane, but a vanishing y component. The mag- 
nitude and sign of the spin polarization are strongly momen- 
tum dependent, and display a complicated interplay between 
singlet-triplet ratio q and SOC strength A. For example, for 
the (,s+p)-wave state, we find that the x-spin polarization in 
the band with the smaller gap can reverse sign near the edge of 
the larger gap. In contrast, the nondegenerate zero-energy flat 
bands of the (d X2y +p)-wave NCS exhibit a strong and robust 
x-spin polarization, whereas the z-spin polarization is smaller 
and changes sign close to ±(kp,- — fcp,+)/2. The doubly de- 
generate states for (c^+^-wave pairing with q > q Cy u give 
opposite contributions to the spin-resolved DOS of unequal 
magnitude, overall leading to a weaker spin polarization than 
for the nondegenerate states. As required by time-reversal 
symmetry, subgap states with opposite edge momenta have 
opposite spin polarizations. In the case of the (s+p)- wave 
state, this ensures the robustness of the linearly dispersing 
edge states against backscattering from nonmagnetic impuri- 
ties. On the other hand, in the (d xy +p)-wave NCS this leads to 
a considerable suppression of time-reversal-symmetric scat- 
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FIG. 3. (Color online) Edge current I y as a function of singlet- triplet 
parameter q for various exchange-field strengths h ex in an NCS with 
an exchange field in the leading edge. Top row: (s+p)-wave NCS 
and exchange field along (a) the x axis and (b) the z axis. Bottom 
row: (d xy +p)-w3Lve NCS and exchange field along (c) the x axis and 
(d) the z axis. 



FIG. 4. (Color online) Edge current I y as a function of exchange- 
field strength h ex for various values of the singlet-triplet parameter q. 
Top row: (s+p)-wave NCS and exchange field along (a) the x axis 
and (b) the z axis. Bottom row: (d xy +p)-wave NCS and exchange 
field along (c) the x axis and (d) the z axis. 



tering processes connecting the oppositely x- spin-polarized 
flat bands for k y > and k y < 0. 

The nontrival spin character of the edge states and the re- 
sulting absence of certain scattering processes could be tested 
experimentally using Fourier-transform scanning tunneling 
spectroscopy. Another possibility is to probe the spin polar- 
ization of the subgap states by bringing the NCS into con- 
tact with a ferromagnetic insulator. It is expected that the 
proximity-induced exchange field H ex will lead to a pertur- 
bative correction to the energy of the spin-polarized edge 
states proportional to X^=i ^ex Pn(^> k y ). Since the flat- 
band states at +k y and — k y are oppositely spin polarized, 
the coupling to the exchange field will shift the energy of 
these edge states in opposite directions. To test this expec- 
tation, we show in Fig. [2] the edge band structure in the 
triplet-dominated phase in the presence of an exchange field 
of magnitude h QX = |H ex | =0.4 acting on the leading edge. 
(The corresponding plot for the singlet-dominated phase is 
given as Fig. [7] in the SM.) Indeed, the edge states of both 
the (s+p)-wave and (d xy +p)-w2ive NCS display a linear shift 
in energy, which is to a good approximation proportional to 
k y ), as long as h QX < A . The doubly de- 
generate flat bands that occur for q > q c jj in the (d xy +p)- 
wave NCS split due to the opposite spin polarization of the 
two degenerate states, see SM. 

Edge Currents. — The energy shifts of the edge states due 
to the coupling to the exchange field can give rise to a sponta- 
neous edge charge current. To evaluate this current, we define 
new electron operators in mixed real- and momentum- space 
c nk y a = ~^T x E/c, c ka e- ifc - n , where n represents the layer 
index and L x denotes the width of the ribbon. With this, the 
layer-dependent current operator along the y axis can be ex- 
pressed as j n = - 1 ^ Y, ky Unly + ffk y ] ' Where N V iS the 



number of k y points in the edge BZ and 

3nXy = _t sin ^2/ (^ C nk y t C nkyt + C \ik y \. C nk y \) > ( 3a ) 

July = + A k V ( C L y i C nkyt + C Lyt C nk y i) ' ^) 

The contribution corresponds to nearest-neighbor hop- 

(2) 

ping, whereas y n J k is due to SOC. We calculate the expec- 
tation value of the edge current at zero temperature from the 
spectrum E\^ y and the wavefunctions \ipi,k y ) of 

L x /2 

^ = -; F EEE Em.<U#*>.0- (4) 

y k y 1=1,2 n=l l,Et<0 

We observe that the current operators Eq. ^ are unaffected 
by the presence of the superconducting gaps or the edge; these 
only enter into Eq. ^ through the eigenstates \^i : k y )- 

In Figs. [3] and [4] we plot the edge current I y as a function 
of singlet- triplet parameter q and exchange-field strength h QX . 
Note that the current reverses sign with the direction of the 
exchange field. For both the (s+p)-wave and (d xy +p)-wave 
cases, we indeed find spontaneous currents flowing along the 
edge of the NCS. However, the currents in the two cases are 
dramatically different. For weak to moderate exchange-field 
strengths, the magnitude of the currents in the (<s+p)-wave 
state is always small compared to the (d xy +p)-wave case. The 
current at the edge of the (s+p)-wave NCS grows linearly 
with the exchange field as long as h QX < A , see Figs.[4ja) 
and[4|b). On the other hand, in the (d xy +p)-wave NCS, the 
current exhibits a remarkable deviation from linear behavior: 
At zero temperature, an infinitesimally small exchange field 
is sufficient to generate a large current in the NCS, as seen 
in Figs.Qc) and|4jd). The currents in both the (s+p)-wave 
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FIG. 5. (Color online) Momentum-resolved edge current I y (k y ) 
for various values of singlet-triplet parameter q. Top row: (s+p)- 
wave NCS and exchange field along (a) the x axis and (b) the z axis. 
Bottom row: (d xy +p)-wave NCS and exchange field along (c) the x 
axis and (d) the z axis. The exchange field is h ex = 0.2. The vertical 
gray lines indicate the projections of the edges of the two spin-orbit- 
split Fermi surfaces. 



and (d xy +p)-wave NCSs depend only weakly on q for both 
Q < Qc,l and q > q Cy u- However, in the (s+p)-wave NCS and 
for an exchange field along the z axis, the current changes sign 
as the pairing state is tuned from triplet dominated to singlet 
dominated. In the (d xy +p)-waYe state, on the other hand, the 
current always shows a clear change in magnitude as one goes 
from q < q c ^ to q > q c ,u> see Figs.[3jc) and[3|d). 

To gain a better understanding of the origin of the currents, 
it is instructive to examine the current I y (k y ) contributed by 
states at +k y and —k y , i.e., the even part of the fc^-summand 
in Eq. In Fig. [5J we show the evolution of I y (k y ) with 
the singlet-triplet parameter q for exchange fields along the 
x and z axes. For the (s+p)-wave pairing state, this re- 
veals that the edge current is predominantly due to states with 
\k y \ < kp,+, whereas in the (d xy +p)-W3LYe NCS the current 
originates mostly from states with kp,+ < \k y \ < kp,-- In- 
terestingly, for an exchange field along the x axis the current 
shows no variation with q in this momentum range and is in- 
dependent of the exchange-field strength (not shown). The k y - 
resolved current for the z-polarized exchange field in Fig.[5jd) 
shows similar characteristics, but also a jump in I y (k y ) at the 
location of the sign change of the z component of the spin 
polarization of the nondegenerate edge states, cf. Fig.[TJf). 

Discussion. — The remarkably strong edge current in the 
(d xy +p)-wave NCS directly follows from the coupling of the 
exchange field to the nondegenerate flat bands of the NCS. 
Inspection of Eq. ([3j shows that the hopping part of the 
momentum-resolved current is proportional to the difference 
in the number of states below the Fermi energy at +k y and 
—k y , whereas the spin-orbit part is proportional to the sum of 
the x-spin polarization at these momenta. Thus, even for an 
infinitesimal field, the opposite energy shift of the oppositely 
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The sign of the edge-state spin polarization along the 
exchange-field direction thus plays a critical role in generating 
the current, and explains the jump in the momentum-resolved 
current for the z -polarized exchange field. We observe that for 
the (gL^ +p)- wave NCS, states with \k y \ £ [&f,+,&f,-] only 
contribute to the current through the SOC coupling, while 
in the (s+p)-wave pairing phase the hopping contribution is 
zero for all k y , since in both cases the number of states be- 
low the Fermi energy at +k y and — k y are always equal. The 
splitting of the doubly degenerate flat bands in the majority- 
singlet (d xy + p)-wave NCS can nevertheless produce a size- 
able current, as the shifting of the negatively (positively) spin- 
polarized states below (above) the Fermi energy generates a 
finite x-spin polarization and hence a current via the spin-orbit 
coupling. This accounts for the change in the current between 
Q < Qc,l and q > q Cy u seen in Figs.|4jc) and (d). 

Experimental detection of the edge currents should be pos- 
sible in spite of the Meissner effect, which implies that screen- 
ing currents exactly compensate the edge currents in a large 
sample. However, whereas the edge-current density decays 
into the bulk on the scale of the coherence length £o> the 
screening only builds up over the scale of the penetration 
depth Xl. For an extreme type-II superconductor, character- 
istic of many NCSs 0, the screening currents will therefore 
be suppressed in a sample of width W with £ <C W <C Xl- 
This argument also holds for an engineered NCS. We note that 
an edge current is also induced by placing the NCS in contact 
with a ferromagnetic metal [ 26 1. 

Summary. — We have shown that the subgap edge states of 
an NCS are typically strongly spin polarized, which we have 
illustrated for two model NCS with distinct topological prop- 
erties, specifically the fully gapped (s+p)-wave and the nodal 
(d xy +p)-wave states. We have demonstrated that the pres- 
ence of an exchange field in the outermost layer generates 
an edge current in the NCS. This is most clearly realized in 
the (d xy +p)-wave NCS, where the nondegenerate flat bands 
acquire a weak dispersion that leads to a strong edge current 
and to a singular dependence on the exchange-field strength. 
The remarkable current characteristics are a novel signature 
of these topologically protected states. 
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Supplemental Figures 

In Fig. [6] we plot the edge DOS and the x and z component of the spin-resolved edge DOS for an (s+p)-wave and a (d xy +p)- 
wave NCS with singlet-dominated pairing state (q = 0.75). This should be compared to the triplet-dominated case shown in 
Fig.[T] Both the continuum and subgap states have a strong spin polarization, which is an odd function of edge momentum k y , 
as required by time-reversal symmetry. The nondegenerate zero-energy flat bands of the (d xy +p)-wave NCS show a particularly 
strong and robust x-spin polarization, with a weaker z-spin polarization that changes sign near ±(kp- — fc j p >+ )/2. The doubly 
degenerate states of the (d xy +p)-wave NCS, which appear in the region \k y \ £ [kF,+,kp,-], give opposite contributions to the 
spin-resolved DOS of unequal magnitude. They therefore exhibit a weaker overall spin polarization than the nondegenerate 
states. For the same value q = 0.75, i.e., also for the singlet-dominated pairing state, Fig.[7]displays the energy- and momentum- 
resolved edge band structure in the presence of a proximity-induced exchange field H ex in the leading edge. The coupling to the 
exchange field shifts the nondegenerate flat bands in opposite directions for k y > kp,+ and k y < — fcp,+ and splits the doubly 
degenerate bands for \k y \ < kp,+- 
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FIG. 6. (a) and (b) Energy- and momentum-resolved DOS 
on a log scale for the ten outermost layers at the (10) edge of 
an (s+p)-wave and (d X y+p)-wayQ NCS, respectively, (c), (d) x 
component and (e), (f) z component of the spin-resolved DOS 
for the outermost layer at the (10) edge in units of h/10 on a 
linear scale. We present typical results for the singlet-dominated 
pairing state with q = 0.75. 



(a) H x lle x s+p (b) H x lle z s+p 



0.2 
tq 0.0 
-0.2 
-0.4 

0.4 
0.2 
tq o.O 
-0.2 
-0.4 



; j 


L J: 





_ ^^^^^ ^^^^ _ 

i i i i i 
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FIG. 7. Energy- and momentum-resolved edge DOS on a log 
scale for the ten outermost layers at the (10) edge of an (s+p)- 
wave and (d xy +p)-wa\e NCS in the presence of an exchange 
field in the leading edge with magnitude h ex = 0.4. (a) and 
(c): exchange field along the x axis, (b) and (d): exchange 
field along the z axis. We present typical results for the singlet- 
dominated pairing state with q = 0.75. 



